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Abstract. We consider a coagulation model first introduced by Redner, Ben- 
Avraham and Krapivsky in [11], the main feature of which is that the reaction 
between a j'-cluster and a fc-cluster results in the creation of a \j — fc|-cluster, 
and not, as in Smoluchowski's model, of a (j'+fc)-cluster. In this paper we prove 
existence and uniqueness of solutions under reasonably general conditions on 
the coagulation coefficients, and we also establish diffcrcnciability properties 
and continuous dependence of solutions. Some interesting invariance proper- 
ties are also proved. Finally, we study the long-time behaviour of solutions, 
and also present a preliminary analysis of their scaling behaviour. 



1. Introduction 

Among the diverse mathematical approaches to modelling the kinetics of clus- 
ter growth, one that has received a good deal of attention consists in the mean 
field models of coagulation-fragmentation type [10] of which Smoluchowski's co- 
agulation system is a prototypical case. The basic dynamic process modelled by 
Smoluchowski's coagulation is the binary reaction between a j-cluster (a cluster 
made up of j identical particles) and a fc-clustcr to produce a (j + /c)-cluster. So, 
the mean cluster size in these coagulating systems tend to increase with time. A 
contrasting case to coagulation is fragmentation in which the basic dynamic process 
is the disintegration of a given j-cluster into two or more clusters of smaller size. 
In these fragmentation systems mean cluster size decreases with time. 

A coagulation system that, in spite of its basic mechanism being binary cluster 
reactions, has cluster size evolution similar to that of a fragmentation system is the 
cluster eating equation. This model was introduced by Redner, Ben-Avraham and 
Kahng in [11] (see also [7]) but has received scant attention since. We shall call it 
the Rcdner-Ben- Avraham-Kahng system (RBK for short) . 

The basic process is the following: when a j-cluster reacts with a fc-cluster, the 
result is the production of a \j — fc|-cluster (see Fig. 1). 

This process is reminiscent of the coagulation-annihilation model with partial 
annihilation in which two or more species of clusters, A and B say, are present, 
and if a cluster Aj reacts with a cluster B% the resulting cluster has size \j — k\ 
and is an A cluster if j > k, is a B cluster if j < k, and is an inert cluster, neither 
A or B, if j = k. Although many problems remain open concerning coagulation 
systems with Smoluchowski's type dynamics and with either complete or incomplete 
annihilation, there is already a relatively large literature about those systems (see, 
for instance, [6, 8, 13] and references therein). Having just one type of clusters, the 
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Figure 1. Schematic reaction in the cluster eating RBK model 



RBK model could result in a somewhat easier system to handle mathematically and 
it is somewhat surprising that, to the best of our knowledge, it has not attracted 
further attention. 

Having in mind the process illustrated in Figure 1, the index j can no longer 
represent the total amount of particles in a cluster (which is a quantity that should 
be conserved in each elementary reaction) but represents only the number of parti- 
cles in a cluster that are, in some sense, active (a concept whose physical meaning 
we must leave undefined.) So, in this paper, every time we refer to a j-cluster, we 
mean a cluster made up with a number j of active particles. 

Assuming the mass action law of chemical kinetics, the rate of change of the 
concentration of j-clusters at time t, Cj(t), has contributions of two different types. 
It decreases due to reactions of the type (j) + (k) — > (j + k), for k G N, which 
corresponds to 

where aj^ is the rate coefficient for these equations. And increases due to 
reactions like (j + k) + (k) — > (j), again with k £ N, which have contributions of 
the type 

^j-\-k,k^j-\-k 

Adding all these contributions we obtain the Redner-Ben-Avraham-Kahng co- 
agulation system 

, oo oo 

(1-1) -ijr = ^2 aj+k,kCj+kCk - ^2 aj^CjCk, j G N. 

fc=i fc=i 

To simplify notation, we shall often write W p . q instead of a Viq c v c q . Naturally, we 
shall always assume that the rate coefficients are symmetric and nonnegative: 

aj,k = a k,j, aj,k > 0, Vj, k e N. 

In this paper we study the existence and uniqueness of solutions (1.1) in appro- 
priate sequence spaces, we investigate some invariance properties of solutions, and 
also start the study of their long-time and scaling behaviours. 



2. Preliminaries 

The mathematical study of (1.1) requires the consideration of appropriate spaces. 
As is usual in works in this area, we will consider the Banach spaces 

oo 

:={x= ( Xj ) G M N | := £i"|ay| < oo }, fi>0 

3 = 1 
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and their nonnegative cones A+ := {x G X^ \ Xj ^ 0}. Observe that the norm 
||c||o of a given cluster distribution c measures the total amount of clusters that are 
present. In the usual coagulation, or coagulation-fragmentation, equations with 
Smoluchowski's coagulation, the norm ||c||i measures the total density or mass of 
the cluster distribution c. Now, with j measuring only the number of active particles 
in a cluster, this norm measures something like an active density or mass. We shall 
omit the word "active" in what follows. 

Being (1.1) an infinite dimensional system, we need some care in defining what 
we mean by a solution. In this paper we use a definition of solution analogous to 
the one in [1] for the standard coagulation-fragmentation: 

Definition 2.1. LetT £ (0,+oo].j4 (mild) solution of the Cauchy problem for (1.1) 
on [0, T) with initial condition c(0) = cq G Xf is a function c = (cj) : [0, T) — > X^( 
such that 

(i) : each Cj : [0, T) — > K. + is continuous and sup ||c(t)||i < oo. 

te[o,T) 
ft °° 

(ii) : for all j G N, and all t G [0, T), we have 



r — . 

/ ^ CLj,kCk(s)ds < OO. 

J ° fe=l 



(iii): for all j G N, the following holds for each t G [0, T), 



c j {t)=c j {0) + 



^2 aj+k,kC j+ k (s)c fe (s) - ^2 a-j.kCj (s)c k (s) 

■fc=l k=l 



ds. 



Remark 2.2. Assuming that, for some nonnegative constant K and all positive 
integers j and k, the rate coefficients satisfy aj^ ^ Kjk, then (i)=>(ii). 
The definition of solution implies that, if c is a solution on [0,T) 7 then each Cj is 
absolutely continuous, so that equation (1.1) is satisfied by c a.e. t G [0,T). 

Also as in works on the standard coagulation-fragmentation equations, we find 
it convenient to consider finite dimensional systems that approximate the infinite 
dimensional equation (1.1). This will be particularly relevant for the existence 
result. 

We will now define the finite dimensional approximation of (1.1) to be considered 
in the paper. In fact, and in contradistinction with the case of Smoluchovski's 
equation, we will prove that, for initial data with compact support, (1.1) reduces 
to this particular finite dimensional system exactly, and so, for that type of initial 
data, the finite dimensional truncation is not an approximation at all but the exact 
system. In fact, it is exactly this compactly supported cases and the corresponding 
finite dimensional systems that Redncr, Ben-Avraham and Kahng considered in 
[11]. 

To motivate the finite dimensional system, consider an initial condition for which 
Cj(0) = if j > N for some positive N. Since the only process is a coagulating one 
in which the resulting cluster has a smaller size, no clusters with size bigger than 
N can be created. Mathematically, this is translated in the finite N- dimensional 
system, for an arbitrarily fixed positive integer N: 

, N-j N 

(2.1) ^ = E WW(<0 - E^.*( c )< 3 G{1,-- -,N}, 

k=l k=l 

where the first sum is defined to be identically zero when N = 1 or if j = N. 
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Naturally, since (2.1) is a finite dimensional system with a polynomial right 
hand side (as a function of the components Cj of the solution vector c = (c? (•))), 
the existence of local solutions to the Cauchy problems follows immediately from 
the standard Picard-Lindclof existence theorem. In the following proposition we 
collect basic results about solutions to this finite dimensional system. 

Proposition 2.3. Let c = (c i (-)) : 7 max -)■ R N be the unique local solution of 

(2.1) with initial condition c(0) = cq and let 7 max be its maximal interval. Then 

(i) : For every sequence (gj), and every m € {1, . . . , N} the following holds 

N 

(2.2) Y. 9i & i = ~ - 9j-k)W jtk - J2 9jWj, k , 

j=m Ti T 2 

where 

T x = T x {m,N) := {(j,k) e {m,...,N} x {1,...,N}\ k^j-m}, 
T 2 = T 2 (m, N) := {(j, k) e {to, . . . , N} x {1, . . . , JV}| k > j - m + 1}. 

(ii) : If all components of the initial condition cq are nonnegative, then also 
Cj(t) ^ 0, for allj e {1,. . .,N} and all t e where R+ := [0, +oo). 

(iii) : sup/ max = +oo. 

Proof. The proofs of these results follow the corresponding ones for the usual 
coagulation-fragmentation equation closely (see [1, Lemmas 2.1, 2.2].) 

(i) : Multiplying equation (2.1) by gj and summing in j from to to N one gets 

N N-lN-j N N 

j—m j—m k—1 j—m k—1 

and now an easy manipulation (a change of notation in the first sum and 
the separation of the second sum into a sum over T\ and another over T-i) 
gives the result (2.2). 

(ii) : Write (2.1) as bj = Rj{c) — Cjifj(c) where 

N-j N 
k=l k=l 

Suppose that, for some r £ I fl K + and all j = 1, . . . , N, we have Cj{r) ^ 
and c r (r) = for some r. For e > consider the initial value problem 



c 



C j( T ) = c j( t )- 

Thus, c%(t) = R r (c e (T)) + e > and, for some rj > 0, c e r {t) > 0, for 
t <E (t,t + t]), for each e G (0, £o), with e fixed. By continuous dependence, 
by making £ j 0, we conclude that c e r (t) c r (t),t 6 [r,r + ??), and thus 
c r (t) ^ 0,t £ [t, r + rj). Hence, we conclude the nonncgativity of each Cj for 

(e/ni+. 

(iii): Using the expression (2.2) proved in (i) and the nonnegativity of solutions 
in (ii) we conclude that, for every nondecreasing positive sequence (gj), we 
have, for all m £ {1, ... , N}, 

N 
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But then c is bounded and, being the right-hand side of (2.1) bounded in 
bounded subsets of R N , we conclude that sup/ ma x = +00. 
This concludes the proof. □ 

It is easy to conclude from the equations (1.1), (2.1) and the Definition 2.1, that 
if c N (•) is a solution of (2.1), then the function c := (c^, c 2 , . . . , cjy, 0, 0, . . .) is a 
solution of (1.1). 

Also easy to conclude, by choosing gj = j in Proposition 2.3-(i), is the fact that 
the density of solutions to the finite dimensional systems (2.1) decreases with time, 
i.e., for all t 2 ^ t\, we have 1 1 c(*2 ) 1 1 1 ^ ll c (^i)l|i- Given the type of coagulation 
process under consideration (remember Figure 1), this is a physically reasonable 
behaviour for solutions to the infinite dimensional system (1.1). However, it is 
not presently clear that other types of (nonphysical) solution can not, in fact, 
exist, similarly to what happens in the case of the pure fragmentation equation [1, 
Example 6.2]. 

We end this section by introducing the following definition: 

Definition 2.4. Let c be a solution to (1.1). 

(i) : we call c an admissible solution if it can be obtained as the uniform limit 
in compact sets o/[0, 00), as N — > 00, of a sequence of solutions c N to (1.1) 
such that Cj 1 = 0, Mj > N. (In particular, (c^ , . . . , cjy) can be a solution to 

^ (2.1).) 

(ii) : we call c a density nonincreasing solution if, for all ti t\, it holds 
||c(t 3 )||i<||c(fi)||i. 

Remark 2.5. In the literature of coagulation-fragmentatiom equations an admissi- 
ble solution is one that can be obtained as the weak limit asn —¥ cc of a sequence to 
finite n-dimensional truncations of the system [3]. In Definition 2.4 we impose the 
condition of uniform convergence in compact subsets of time t. This corresponds 
to what we can prove in this case (see Corollary 3.2 below); it is also what hap- 
pens to be the case for coagulation-fragmentation equations with coagulation kernels 
growing at most linearly [1, Corollary 2.6]. 

3. Existence of solutions 

In this section we shall prove existence of solutions in X* of Cauchy problems 
for (1.1) with initial data in X± , with some mild conditions on the coagulation 
coefficients. 

Theorem 3.1. Assume a,j t k ^ Kjk, for some positive constant K and all positive 
integers j and k. Let cq G X^ . Then, there is at least one solution of (1.1) with 
initial condition c(0) = Co, defined on [0, T), for some T S (0, +00]. 

Proof. As is usual in coagulation studies [1, 9] the proof is based on passing to 
the limit N — > 00 in a sequence of solutions to the iV-dimcnsional system (2.1), 
which we do by an application of Helly's selection theorem, and then by proving 
that the limit sequence is a solution to the infinite dimensional system (1.1). In 
order to do this we need some bounds on the moments of the solutions to the finite 
dimensional systems. Actually, for the RBK system, the application of the method 
just described is much easier than in [1, 9] for the coagulation- fragmentation system 
with Smoluchowski coagulation due to the a priori estimate (3.1) on the zeroth and 
first moments. 
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Let c N be the solution to the finite iV-dimensional system (2.1) satisfying the 



initial condition (0) = cqj, for j 6 {1, . . . , N}. By putting gj = 1 and gj = j 
(2.2) we immediately conclude that, for 1 ^ m ^ N, 



in 



JV 



N 



£ef<0, and X • °> 
respectively. Thus, for p = 0, 1 we have 

JV JV JV oo oo 

(3.i) Y,f c ? < Y,i p 4 = E^ < E^ < E^ = imi 



I?- 



j=m 



j=m 



j=rn 



3=1 



Let us now prove that (c w ) is uniformly bounded in W 1 ' 1 ^,! 1 ), for all fixed T E 
(0, oo). From the definition of the norms in Xq and X\, and from (3.1) with p = 
we immediately get 

(3 - 2) 

\\cf \\ lH o,t) = [ \cf(s)\ds^ f \\c N (s) \\ ds < / ||c ||o<is = T||c ||o<T||c t ,||i 
Jo Jo Jo 

By equation (2.1) we have 



dc 



.JV 



L!(0,T) 



dt 



» 



(3.3) 



^ TV — j jV 

< / y2w j+k . k (c N (s))ds+ [ Y^W 3 . k {c N { S ))ds 

J ° k=l J ° fe=l 



Estimating the first integral in (3.3) we obtain 



rp j\ — J rp jy — ^ 

/ "£W j+k , k (c N (s))ds < K[ "£jkcy +k (s)c»(s)ds 
J ° k=i Jo k=i 



T N-j N-j 

Ja k=i k=i 

^ K [ \\cP(8)\\ld8 < K [ \\co\Hds 
Jo Jo 

= KT\\co\\l 
and for the second integral in (3.3) we get 

T n T JV 

/ J2 w ^( cN ( s » ds < K 
J ° k=i J ° k=i 

-T JV JV 

< K E^ c f( s )E fcc ^( s ) ds 
J ° k=i 



k=l 



< K ( T \\c N {s)\\\ds < K f T \\co\\\ds 
Jo Jo 

= KT\\c4\. 
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Thus, substituting in (3.3) we conclude that 



(3.4) 

and therefore 



def 



dt 



< 2KT\\co\\l 



L!(0,T) 



IKIUw) + 



dc 



N 



dt 



s=:(l + 2X|| C o||i)T|| C o|h 



So, by Helly's selection theorem, for each fixed j there exists a subsequence of 
(c^)jv (not relabeled), converging pointwise to a BV function in [0,T], Cj(-), 

cf(t) -> c,-(i), as TV -> oo, V* £ [0, T], Vj £ N. 
But then, for each g £ N, and for each t £ [0, T], 

^2jcf(t) -> X] j ' c j(*)' as ^ ^ 00 - 

and therefore, by (3.1), for any such q, Ylj=i3 c j(t) ^ ll c o||i- By making q — > 00, 
we obtain 



(3.5) 



Since proposition 2.3(h) implies Cj(t) ^? 0, this proves that, not only c(t) £ X± , for 
each t £ [0, T], but also that condition (i) of definition 2.1 is fulfilled. 

It remains to be proven that the limit functions Cj solve the RBK system (1.1). 
In order to obtain this result, we shall pass to the limit N — > 00 in the equation for 



(*) = CQj + 



f J2 W J+k , k (c N (s))ds - fY, WjM^Wda. 
Jo k=i Jo k=i 



Thus, we need to prove that, for all t £ [0, T], 



(3.6) 

and 

(3.7) 



t N 
k=l 



N- 



t 00 



J2WjAc(s))ds, 



rt N ~3 



k=l 



The proofs of (3.6) and (3.7) are entirely analogous, and so we shall present only 
the proof of (3.6), leaving the details of the other to the reader. 

We first start by proving that the right-hand side of (3.6) is well defined. Let p 
be an arbitrarily fixed positive integer. By the definition of (cj) we know that 



EN N 
a l,kCj C k 



fc=l 



E 

k=l 



8 



F.P. DA COSTA, J.T. PINTO, AND R. SASPORTES 



and from (3.1) we have that, for all positive integers TV and p, 

f>, fe cfcf <K\\co\\l 
fc=i 

and thus also 

p 

J~] aj,kCjCk < A'||c ||?. 
fc=i 

Consequently, since the right-hand side is independent of p and all the terms are 
nonnegativc, 

oo 

y~]aj,kCjCk < A||c ||i, 
fc=i 

and the dominated convergence theorem implies that, for all t £ (0, T), with T < oo, 
the right-hand side of (3.6) is well defined. 

Now we shall prove the limit in (3.6) holds. Let m be a positive integer such 
that 1 ^ m < N < oo but otherwise arbitrarily fixed. Then 



,t n 



r _ _ r _ _ 

/ Y] a j,kC? (s)c%(s)ds - / yV-, fc Cj-(s)c fc (s)ds 
- 70 fc=i - 70 fc=i 

„ t m— 1 

(3.8) < / Y / ^,k\4(s)c%(s)-c 1 ( S )c k ( S )\d S + 

J ° k=i 



(3.9) 



+ / Y a i' kC j'^ c k( s ) ds + Y aj,kCj(s)c k (s)ds, 

J ° k= m J ° k=m 



and we need to prove that the right-hand side of this inequality can be made 
arbitrarily small when N — > oo, by choosing m sufficiently large. 

Since each term in the sum in (3.8) converges pointwise to zero, the sum has a 
finite fixed number of terms, and its absolute value is bounded above by 2.ftr||co||i, 
the dominated convergence theorem implies that (3.8) converges to zero as N — > oo. 

Let us now consider the integrals in (3.9). Define p m := ||co[|i J2 < jL m j c oj- Clearly 
p m — >• as m — > oo. 

From (3.1) we conclude that 



f £ %)fe cf ( S ) C f ( S )d S ^ K[ £j<$(8)k<$(8)d8 

Jo , Jo , 

k—m k—m 

K ||co||i V kc%(s)ds 4 K Pm ds 
Jo , „ Jo 



k 

(3.10) = KTp m , 

and so we get the first integral in (3.9) can be made arbitrarily small by choosing 
m sufficiently large. For the second integral the result is proved in an analogous 
way: For all 1 ^ m < p we have 



p 

EN N \ " 



fc= 
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Due to (3.1), the sum in the left-hand side is bounded by Kp m , and so we also get 

p 

E a>j,kCjCk < Kp m , 

k—m 

for all p. Since this bound is uniform in p, we have 



E aj,kcjc k — 



p— >oc 

k—m k—m 

Hence, by the dominated convergence theorem, the second integral in (3.9) can also 
be made arbitrarily small by choosing m and N sufficiently large. 

This completes the proof of (3.6). As pointed out above, the proof of (3.7) is 
entirely analogous and will be omitted. □ 

Corollary 3.2. The solution obtained in Theorem 3.1 can be extended to t £ 
[0,+oo[ as an admissible solution. 

Proof. The uniform convergence property is again a consequence of (3.1). In fact, 
by applying to (3.1) an argument similar to the one that led us to (3.5) we obtain, 
for each m, N G N, t G [0,T], 

OO OO 

j—m j—m 

Since X^m-? c oj — ^ 0; as m — > oo, we conclude that the series in the l.h.s. of this 
inequality with m = 1 is convergent uniformly in (i, N) G [0, T] x N. Since, for each 
j and t, j\cj(t) — \t)\ — >• 0, as N — > oo, wc conclude that, as iV — > oo, 

oo oo 
0=1 3=1 

as W — ^ oo, uniformly in i G [0, T]. 

That c(-) is extendable to [0, +oo[ is a consequence of the arbitrariness of T > 
and estimate (3.5). □ 

4. The moments' equation 

As in the studies of the usual coagulation- fragmentation systems, a weak formu- 
lation of (1.1) is a tool of the utmost importance. This weak version, presented 
next, is the version of the expression (i) of Proposition 2.3, written for (1.1) instead 
of (2.1). 

Proposition 4.1. Let c = (<%•(•)) : I max —> M. N be a solution of (1.1) and let 
T,t G /max be such that r ^ t. For every sequence (gj), and all positive integers m 
and n with m < n the following moment's equation holds 



(4- 1 ) E ^ r > U] E •''''•• :r! = 

j=m 

I J2^o ~93-k)w j>k - f E '' 11 -* • / E^-fc^fc- 

Jt a. J t q Jt c_ 



j=m j=m 



S 2 S3 
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whe 



51 = Si(m, n) 

52 = Sz(m, n) 

53 = S 3 (m,n) 



= {(j, k) e N 2 | m + 1 < j < n + 1, 1 < k < j - m}, 
= {(j, k) e N 2 \ m < j < n + 1, k > j - m + 1} 
= {(j, k) £ N 2 \ j > n + 1, j - n < k < j - m}. 



In Fig. 2 we give a geometric representation of the regions Sj 



1- 




FiGURE 2. Regions Sj defined in Proposition 4.1 



Lemma 4.2. Suppose the coefficients aj_k satisfy the condition Qj k ^ Kjk, then 



(4.2) 



lim / m > W^- kds = 0, 



S2("1,CXd) 

Proof. Considering <?j = 1 in the moments' equation (4.1) we obtain 

(4.3) Yl c i (*) - E c *w = - / E + / E 

j— m i=m T S2 T S3 

We start by estimating the expression in the S3 region. We clearly have 

00 j—m 

Y, Wj, k < K Y E Ucj)(kc k ) < AlcWH?, 

S3(m,n) j=n+l k=j — n 

00 

and also \J Wj,fc ^ i£"||c(t)||i jcj(t) —} 0, pointwise as n — > 00. Since, 

S 3 (m,n) j=n+l 

from the definition of solution, ||c(f)||i is bounded in [t, t], applying the dominated 
convergence theorem gives 

/ Y. Wj^ds — > 0, as n —> 00. 

S 3 (m,n) 
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In S2 the same bound Wj,k < K\\c(t)\\ x is true and the convergence 

S2(m,n) 



n 00 



00 00 



E E » -E E ^ 

j=m k=j— m+1 j— m k=j— 

is valid pointwise in t as n -> 00. Hence, again by the dominated convergence 
theorem, 

/ E / E ^' fc ' as n 00 ■ 

J T c f J T a f \ 



S2(m,n) 



Now taking limits, as n -> oo on both sides of (4.3) we obtain 



OO OO 

E^(*)-E c iW = -/ E 



W 



S2(m,oo) 



But ^Y^jLm c 3 ^ S^lmJ c j — >■ as m — > oo, since by definition of solution, 
c G Xi, we then obtain 

(4.4) 

This concludes the proof. 



lim 

m— ^oo 



J m ^ W jlk ds = Q. 



S2{m,oo) 



□ 



Another estimate that will be useful is the following 



Proposition 4.3. Suppose the coefficients a,j,k satisfy the condition dj^ $J Kjk, 

and the sequence (gj) satisfies \gj \ ^ j then, for each m <EN, 

(4.5) 

OO OO „t r 

E SjCjit) ~ E 9]Cj{T) = - lirn^ / ^ - g,-_fc)W,-, fc + 9jWj, k ■ 

j=m j=m r Si(m.n) S2(m,n) 

Furthermore, with the stronger assumptions aj,k *s K{jk)^ with ^ (3 ^ ^, and 
the sequence (<?j) satisfying \gj\ ^ j and \gj — gu\ ^ M\j — k\, for all j and k, and 
for some positive constant M , the following holds true: 

oo oo r-t ft 

(4-6) E 9jCj(t)-J2 9jCj(T) = -J Y iSi-9i-h)W jtk - J E 9jW jtk . 

j=m j=m T Si(m,oo) T S2(m,oo) 

Proof. In the moments' equation (4.1) we prove that 



(4.7) 



J E .'/.. / " ; / — > as moo 



S3(m,n 



In fact, we observe that S^{m^n) C Szin + l,oo), and since in 63(771, n) it holds 
that \j — k\ = j — k n, we get, by the previous lemma, 



0^ 



f E 95-kW^k f X b' ~ fc l^fc " 

&3{m,n) 

^ (n + 1) y 51 ^> fc ^ °' 



S3(m,n) 



as 71 — r OO , 



S2(n+l,oo) 
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thus proving (4.7). As a consequence, by taking the limit as n —> oo in (4.1), we 
obtain (4.5). Now, by imposing the stronger conditions of the second part of the 
proposition, we have, for each n G N, 

n 

(4.8) J2 9jW jtk K E j(jkfc jCk ^K £ jkv.c, ■ K c , ^\,V,. 

S 2 (m,n) S 2 {m.n) S 2 (m,n) j=m 

where in the second inequality we have used the fact that, if (j, k) € S2{m, n), then 
j ^ k, and so {jk) 13 < k 213 < fc, due to the assumption (3 ^ |. Similarly, for (j, fc) 
in Si(m, n) we have k < j, and thus 

(4.9) ^ //, .</., , Wa ■ MK £ Wkfcjck 

S± (m.n) Si(m,n) 

n 

< mk £ i kc i Ck < MK \\ c h £ fa- 

Si(m,n) j=m 

Estimates (4.8) and (4.9) and (i) from the Definition 2.1, together with the domi- 
nated convergence theorem allow us to prove (4.6). □ 

Corollary 4.4. Suppose the coefficients a,j t k satisfy the condition aj t k ^ Kjk, 
then any solution is a density nonincreasing solution. Moreover, with the stronger 
assumption that aj.k ^ K(jk)P with 0^/3^ h, the following holds true 

oo oo 

(4.10) > >,(*) - > >,(r) 2 / X Ma - I > ..J".,,- 



-£**w = - 2 / £ £ ^ - / £j 

3=1 j = l Jr j=lfc=j+l Jt 3 = 1 

OO CXD „f OO OO 1 OO 

(4.1D E^)-E c ^) = 4/ EE^- / E^ 



3=1 3 = 1 3 = 1 fe=l 3=1 

Equations (4.10) and (4.11) follow from (4.6) considering gj = j and gj = 1. 
An interesting particular case concerns the evolution of the number of clusters 
of odd size, that we will consider in Section 9. 

Corollary 4.5. Choosing gj = Sj i0 dd then 

OO OO oo oo 

(4.i2) Yl - E c ^ = / E E 

3 = 1 3 = 1 ^ T 3 = 1 fe=l 

3 odd 3 odd j odd fc odd 

Proof. Considering this choice of the sequence (gj) in (4.6) we have, after some 
rearrangements 



£<*(*)- $>j(r) = -/ £^-/ £%*+/ E 1 

3 = 1 J = l Jt S!US 2 Jt S 2 Jt Si 

j odd j odd j odd j odd j even 

k odd fe even k odd 



3,fc" 



The last two terms cancel out, since using the fact that Wj t k = Wk,j we have 

oo oo oo oo oo oo 

E = E E = E E ^ = E E ^ 

S 2 3=2 fc=l fc=l 3 = fc+l 3 = 1 fc=3+l 

3 odd 3 odd k even fe even j dd 3 even fe odd 

k even 
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and similarly in the Si region 



^ Wj,k = V] = Wj ik , since j is even and fc is odd. 

Si 3=1 fc=j i=i fc=i+i 

j even j even fc odd j even odd 

A: odd 

oo oo 

Since Wj.fc = /J Wj'.fcj this concludes the proof. □ 

S lU S 2 j=l fe=l 

j odd j odd fc odd 

fc odd 

5. A UNIQUENESS RESULT 

Wc now consider a uniqueness result for (1.1). The result is obtained by assuming 
the initial value problem has two solutions and proving they are equal. This will be 
done, as usual in coagulation problems (see, e.g., [1, 2]) by appropriate estimates 
on the solutions and the use of Gronwall's inequality. The proof requires conditions 
on the coagulation coefficients that are slightly more restrictive than the ones used 
for the existence result. At the time of writing it is not clear if these conditions can 
be significantly relaxed. 

Proposition 5.1. Let ^ K(jk)^ , with (3 ^ i. Then, for each c <E there 
is one and only one density nonincreasing solution in [0,T) such that c(0) = cq. 

Proof. Suppose the initial value problem for (1.1) with the initial condition c(0) = 
Co G X^ has two density nonincreasing solutions, c and d. Let x(t) := c(t) — d(t) 
and Mj t k '■= aj^icjCk — djdk) = a,j_k(c.jXk + dkXj). We shall prove that c = d by 
establishing that, for some a, the sum j a \xj\ is identically zero. This will 

be achieved by deriving an inequality for this quantity and applying Gronwall's 
inequality 

So, let us consider Proposition 4.1 with m = 1. From the above definitions we 
get 

X! 9i x i = I L ~ 9j) M 3,k ~ X •'/./ A + 3j-kM^ k j ds. 

J = l J ° V Si s 2 s 3 ' 

For each t 6 [0, T) consider gj ~ j a sgn(ccj) for some a ^ j3 such that a + (3 ^ 1. 
We now estimate the sums over each of the Sj. For (j, k) <E S\ we get 

(9j-k - 9j)x k = ((j - k) a sga(xj- k ) - j a sgn(xj))x k 

= {(j - k) a sgn(xj- k ) - j Q sgn(xj)) sgn(a; fc )|a; fc 
= ((.? - k T sgn(xj- k x k ) - j a sgn(x jX k )) \x k \ 

< (U - k) a + 3 a )\x k \ 

< 2j a \x k \ t 

and, by a similar computation, (gj_k~9j)%j = k) a sgn.(xj-k)—j a s S u ( x j)) x j — 
(0 — k) a sgn(xj-k%j) — j a ) \xj | < 0. Using these bounds and the assumptions on 
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a and /3 we can estimate the sum over Si as 

n t n j—l 



J ° Si J ° J=2fe=l 

< 2#||c ||i / YV|a*|tfc. 



For (j, k) € 5*2 we have — j a sgn(xj)x k ^ j Q |£A;|, and j Q sgn(a;j)xj = j Q |a;j|, from 
which it follows that 

- J29oM jtk ds < K Y,Y.3 a ^fi.^k\-d k \ Xj \)ds 

J ° S 2 J ° j=lk=j 

ft oo 

^ A'Hcoiii / yvisfcids. 
J ° k =i 

Finally, for (j,k) € S3 we have the estimates (j — k) a sgn(xj- k )xk ^ n a \x k \, and 
(j — k) a sgn(x j)xj = n a \xj\, from which it follows 

ft j-t 00 j-i 

J29j-kMj, k ds ^ Kn a (jk) p {cj\x k \ + d k \ Xj \)ds 

J ° S 3 J ° ]=n+l k =j-n 

(5.1) < Kn a I 3 Pc i E ^\x k \ds + 

j—n+l k —j — n 

(5.2) +Kn a J2 3 P \ x i\ E k0d kds- 

J° j=n+l k=j-n 

For the double sum in (5.1), using a + f3 $J 1 we conclude that 

j—l 1 00 j—l n I, 00 

, a, , , II c II x , «, 



OO J — 1 _.00 j—l j, II 00 

5: f Cl E E E fcV-K^E^ 

j'=n+l k—j—n j—n+l k —j — n k—1 

For the double sum in (5.2) we have, again using a + /3 $J 1, 

OO j—l || || / OO OO s 

E /N E E A* + E /<H 

j— n+l k=j—n ^ j—n+l j—n+l 

and therefore 

/ E 5j _ fc M Jifc( is < X||c ||i / ^k f} \x k \ds + 
Jo S 3 70 k=i 

+^iicoiii / E A* + E i" d iK 

Combining the estimates on the three regions we get 
(5-3) 

]TrN< / ^/K-i + ^Hcolli / ( E + £ fd, 

j=l Jo j=l Jo \=n+l J=n+1 
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By the definition of solution and the assumption ft < 1/2 < 1, 

oo oo 

J^Cj ^ jcj — > 0, as n — >• oo, pointwisely, 

and since c(-) is a density nonincreasing solution, 

oo oo oo 

J=n+1 j=n+l 3=1 



Thus, by the dominated convergence theorem 

/ j^Cjds — >• 0, as n — ► oo, 

^° j=n+l 



the same being valid for Y^jL n +i J^dj. 

Therefore, letting n — > oo in (5.3) and using the assumption a ^ ft in the 
right-hand side we obtain 



oo „t oo 

Ei a N<4iq|co||i / ]T?>^ S 

Jo 



3=1 " u 3=1 

Since Xj(0) = 0, by Gronwall inequality we conclude 

oo 

Ei a K-| = o, We[o,T), 

3=1 

and so Xj = 0, for all t £ [0, T) and j £ N, thus proving uniqueness. □ 

Remark 5.2. Observe that, if ft > 1/2, it is not possible to find numbers a such 
that a + ft ^ 1 and a ^ ft. It is this elementary technical reason that forced us 
to consider ft ^ 1/2 in the uniqueness result in Proposition 5.1, since in this case 
such numbers a obviously exist. It is not presently clear if a uniqueness result is 
true in more general situations. 



6. Differentiability and continuous dependence 

In Theorem 3.1 and Corollary 3.2 we proved that with the hypothesis that a,j k ^ 
K jk, for each initial condition Co £ X± , there exists at least one admissible solution. 
With the stronger assumption that a,^ ^ K(jk) 1 ^ 2 , Theorem 5.1 implies that there 
is a unique solution defined in [0, oo) and so, it has to be an admissible solution. 
We now address the issue of differentiability of such solutions. 

Theorem 6.1. If there is K > such that, for all j,k £ N, aj,u ^ Kjk, and 
c(-) = (cj(')) is an admissible solution then, the functions t <— > Cj(t), for j £ N, and 
t i — ^ Sj°=i c j{t) are continuously differentiable. Moreover, (1.1) is satisfied for all 
t £ [0, oo) and 

d ' ^ 

(°) iiE^w = "EEwW ct W' vi[o,+oo). 
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Proof. If a solution is admissible, we can use (3.1), like in Theorem 3.1 and Corol- 
lary 3.2, to conclude that Yl'jLi J c o ft) ^ s uniformly convergent in compact subsets 
of [0, oo). This, together with the assumption on the coefficients aj ik and the con- 
tinuity of Cj(-), for j G N, allows us to prove the continuity of the right-hand side 
of (1.1) in t, thus proving the continuous differentiability of each c 3 (-), and the fact 
that (1.1) is satisfied by c(t), for all t G [0, +oo). 

From the proof of Lemma 4.2, we already know that 

oo oo p f oo k 

E c o ft ) ~ E c o ( T ) = - EE a i> kC i ft) Ck ( s ) ds ■ 

j=l j=l Jr fc=l J=l 

Since the double series of continuous functions in the r.h.s is uniformly convergent 
in each compact interval, the result follows. □ 

With an extra condition on the kinetic coefficients we can say more about the 
differentiability of the moments for any solution: 

Proposition 6.2. If there is K > such that, for all j, k G N, a,j t k ^ K{jk)P , with 
(3 ^ i, and c(-) = (cj(-)) is a solution of (1.1), and (gj) is a nonegative sequence 
such that ^ gj ^ j, and, for some positive constant M, \gj — gj-k\ ^ Mk for 
1 ^ k *S j, then for m € N, i 4 Tl C jLrn9j c ]ft) * s continuously differ entiable and 
moreover, for t 0, 

d °° 

( 6 - 2 ) j t J2^ft) = - E fa - u-MMt)) - E 9oWMt)). 

j—rn Si(m,oo) S2(m,oo) 

In particular the differential versions of (4- 10), (4-H) cmd (4-12) hold true. 

Proof. From the estimates (4.8) and (4.9) taken together with the property of 
uniform convergence of Y^jLi3 c 3> ^ i s clear that the convergence of the double 
series in (4.6) is uniform so that the series in the integrals of the r.h.s. of that 
expression define continuous functions in t. This implies our result. □ 



With respect to the continuous dependence relatively to the initial conditions 
we prove the following partial result: 

Proposition 6.3. If < K(jk)P wth f3 ^ |, if a + (3 ^ 1 and if c and d are 

solutions of (1.1) satisfying c(0) = Co and d(0) = do then, for each t ^ 0, there is 
a positive C(t, ||co||i) such that 

(6.3) \\c(t) - d(t)\\ a sc C(t, ||co||i)||c - do|U. 

Proof. By writing 



Cjft) = c 0j + 



dj(t) = doj + 



22a j+k ,kCj+k(s)ck(s) -^2a jtk Cj(s)c k (s) 

k=l k=l 

oo oo 

y] aj + k,kdj + k(s)d k (s) - y^aj,fcdj(s)d fc (s) 



ds, 



ds. 



k=l 



k=l 
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and defining x(t) = c(t) — d(t) we perform the same estimates as in the proof of 
Proposition 5.1 to obtain, this time, 



-J CO WW 

^iicoiii / E A J + xiMoiii / E A 



3=1 3=1 " u 3=1 

-t oo »t oo 



3 ' 

3=n+l " u 3=n+l 

instead of (5.3). Hence, by making n — >■ oo and using the same arguments as in 
that proof, we obtain, 

oo oo r i oo 

"Iff,- 



Ej' q i^wkEj' q i^(°)I + 4A 1 c oIIi / E^' Q i 

3=1 3=1 3=1 

By using the Gronwall lemma estimate (6.3) follows. □ 



Remark 6.4. Here we recall the Remark 5.2 with respect to our restricting hypoth- 
esis on the kinetic coefficents and on the growing rate of (gj) . Notice that we were 
not able to prove continuous dependence on the initial conditions with respect to 
the X\ norm, except in the case where a,j,k K, for all j, k £ N. At the moment it 
is not clear to us whether this is an essential feature of the RBK equation or if it 
is just a technical limitation due to the methods we have used. It is possible that, 
in a more general case, instead of a continuity property based on a norm estimate 
like (6.3), we can prove an upper semicontinuity property for admissible solutions 
similar to that in [1, Theorem 5.4]. 

7. Some invariance properties of solutions 

From the physical process that we arc modelling (sec scheme in Fig. 1) it is 
natural to expect that if initially there are no clusters of size larger than p, then 
none will be produced afterward. This is established next. 

Proposition 7.1. Assume the Cauchy problems for (1.1) have unique solutions. 
Then, for every pgN, the sets 

:={ceX+\c j = 0,Vj>p} 

are positively invariant for (1.1). 

Proof. Let c be a solution of (1.1) such that c(t) = Co € X^ p , for some r ^ 0. 
Let c p () be the unique solution of the p-dimcnsional Cauchy problem 

EH--:,-"; £V M (c") 

fe=l fe=l 

for j = 1, . . . ,p (with the first sum defined to be zero if j = p.) Then, the function 
(c\ , cf, . . . , c p , 0, 0, . . .) is a solution of the infinite dimensional system (1.1) and, by 
uniqueness, it must be the solution c. Therefore, for all t ^ r, we have Cj(t) = 
when j = p + l,p + 2, . . ., that is, c(t) £ X^ p for all t ^ r, which proves the 
result. □ 
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This invariance property also occur in fragmentation equations: if the initial 
distribution of clusters has no clusters with size larger than p, then they cannot 
be produced by fragmentation of those (smaller) ones that are initially present; 
a reasonable enough result. Invariant sets for coagulation equations with Smolu- 
chowski coagulation processes arc of a different kind but one can also characterize 
them without much difficulty [4], In fact, we can use a similar approach also in this 
case to characterize the positivity properties of the cluster distribution (i.e., the 
subscripts j for which Cj(t) > 0) in terms of those same properties for the initial 
data. Let us first introduce some notation. 

For a solution c = (c,-) to (1.1), denote by P := {j G N|c,(0) > 0} the set of 
integers (finite or infinite) describing the positive components of the initial condition 
(cj (0)) , and let gcd(P) be the greatest common divisor of the elements of P. Define 
J{t) := {j G N|cj(t) > 0}, the set of indices for which the component of the 
solution is positive at the instant t. Naturally, P = J'(O). Now we have the following 
result: 

Proposition 7.2. Assume uniqueness of solution to initial value problems for (1.1) 
holds. Let #P > 1. Then, 

gcd(P) = m => J(t) = mNn [l.supP], Vi > 0. 

Proof. We first remark that, by uniqueness of solution and the form of the system 
(1.1), if, for any j, one has Cj(s) = for all s in a nondegenerate interval [r, t], then 
Cj is identically zero for all times. 

In order to prove the proposition it is convenient to write system (1.1) in a dif- 
ferent form, similar to what was done for the Smoluchowski's coagulation equation 
in [4]. Define 

oo oo / /"* \ 

RjW) ~ y^, a j+k,kCj+k(t)ck(t), (Pj(t) := y^fcCfcft), Ej(t) := exp / (pj(s)ds). 

fe=i fc=i VJo ' 

Observe that all these functions are nonnegative, for nonnegative solutions and for 
all t (and furthermore Ej(i) ^ 1). Using them, write (1.1) as 

Cj = Rj - cm 

and apply the variation of constants formula to get, for all i r ^ 0, 
(7.1) Cj(t)Ej(t) = c 3 {t)E 3 {t) + J E j (s)R j (s)ds. 

Equation (7.1) allows the following conclusions to be immediately drawn: 

(i) : if Cj(t) > 0, then Cj(t) > 0, for all t > r, or, equivalently, J{t) D J{t) 
for all t > r, and in particular J(t) D P, for all t > 0. 

(ii) : due to the definition of Rj(-), if li > t\ are two numbers in J(t), then 
li-t\ € J(t), for all t > t. 

(iii) : by (i) and (ii) one concludes that, if p%, . . . ,p n G <J(t), (assuming, with- 
out loss of generality, that pi < pj for i < j), then, for all integers m\ G Z, 
we have m\p\ + . . . + m n p n G Jit\ for t > r, provided the integers Wj are 
such that 1 ^ rriipi + . . . + m n p n ^ p n . 

Suppose that 1 < #P < oo. Let us write P = {p%, . , . ,p n }. By Bezout's lemma 
in elementary number theory [12, Chapter 1] we conclude that if gcd(P) = m 
then the smallest positive value of mipi + . . . + m n p n is m and all other larger 
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values are multiples of m. Hence, this result, (iii), and Proposition 7.1 imply that 
J(t) D mN n [1, supP], for all t > 0, when the initial data is finitely supported. It 
is clear that, if m = 1, then equality holds and the proof is complete for finitely 
supported initial data. 

In order to complete the proof we now need to prove that, when m > 1, we also 
have J(t) C mN n [1, sup P] . 

Wc first note that gcd(P) = m =>• P C mN (~l [1, sup P] . This is obvious since the 
assumption implies that every element of P is a multiple of m (hence an element 
of mN) and, naturally, it is not bigger than the supremum of P. 

Let us now prove the result. Note that it suffices to prove that, for any q 6 N, 

q mN n [I, sup P] => c q (t) = 0, Vt > 0. 

Let d = (dj), j 6 N, be defined by dj = if j > supP, and, for j ^ supP, let dj be 
given by the solution of the ordinary differential equation 

{dj = if j gmNn [l,supP], 

p-j oo 
dj = a j+k.kd J+ kdk - d 3 a j,kdk otherwise, 
fe=l k=l 

where p := supP, with initial condition 

d,-(0)=c J -(0), ,j = l,...,p. 

Let j <j£ mN n [l,supP]. We know that j P and therefore 4, (0) = Cj(0) = 0. 
Thus, by the differential equation, dj(t) = 0, for all t ^ 0. On the other hand, if 
j ^ mNn [1, sup P] it is not possible that both k and j + k belong to mND [1, supP], 
for every k. Thus, we conclude that 

p-j 

/] Q>j+k,kdj+kdk = 0. 
fc=i 

Moreover, since dj = for all t, we also have 

p 

d j^2 a j-k d k = 0, 

fc=i 

and hence d is also solution of the system. 

p-j p 
dj = aj + k,kdj+kdk ~ dj ^aj ;fc d fc , j = 1, . . . ,p, 

k=l k=l 

with initial condition dj(0) = Cj(0). Therefore, by uniqueness, c = d and we con- 
clude that q mN n [1, supP] c,(t) = d 9 (t) = 0, Vt^ 0. 

Suppose now that #P = cxd. Let gcd(P) = m. By [5, Proposition 5] there exists 
a finite subset P n C P, with #P„ > 1, such that gcd(P„) = m. Apply now the 
above argument to P n instead of P. This concludes the proof. □ 

The previous proof does not hold if #P = I. In that case, a peculiar behaviour 
occurs, not exhibited either by the usual Smoluchowski's equations, or by the frag- 
mentation equations. 
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Proposition 7.3. Assume uniqueness of solutions to initial value problems for 
(1.1) holds. If a solution to (1.1) starts monodisperse, it stays monodisperse for all 
later times, or, in the notation used above, 

#P = 1 => J(t) = P, Vi > 0. 

Proof. Let Cj(0) = X5j iP , A > 0, for some positive integer p. Then, c(0) £ X^ p and, 
by Proposition 7.1, c(t) £ X^ p , for all t > 0. Therefore c 3 (t) = for all j ^ p + 1 
and, for j = 1, . . . ,p, Cj(i) is given by the p-dimensional system considered above. 
Now let Cj(t) = a(t)8j^p. Obviously we have, for all j = 1, . . . ,p, 

p-j p 

W 3+k . k {c) = 0, W ^(c) = a j>p a{t) 2 S jiP . 

k=l k=l 

Thus, c(-) will be a solution of (1.1) with initial condition {X8j, p ) if and only if a(-) 
solves 

{a = —ap :P a 2 
a(0) = A. 

Hence, solving this initial value problem and substituting back into the expres- 
sion for c(-) we obtain the following solution of (1.1) 

= m — 7 

By uniqueness, it is the only solution satisfying the initial condition Cj(0) = XSj tP , 
which proves the result. □ 



8. ON THE LONG-TIME BEHAVIOUR OF SOLUTIONS 

In this section we start the investigation of the long-time behaviour of solutions. 

Having present the physical process under consideration (see scheme in Fig. 1), 
it is natural to expect that, under rather mild conditions, all solutions will converge 
pointwise to zero as t — > oo. For solutions in X\ this can be phrased by saying that 
solutions converge to zero in the weak-* sense [2, page 672]. 

The result is given in the following proposition, the proof of which is rather 
easy and follows the same ideas used to prove the same result in the fragmentation 
equation [3, Theorem 4.1], which we reproduce here for the sake of completeness of 
presentation. 

Proposition 8.1. Let dj t k ^ Kjk, and let c be a solution of (1.1) with initial 
condition c(to) = Co £ X^~ , defined on [io,oo). Assume that ajj > 0, for all j. 
Then, for all j g N it holds that Cj(t) — > as t —> oo. 

Proof. Let t > r, and consider the moments' equation (4.5) with gj = 1, 

oo oo „i n oo 

E Cj(*) - E c j( t ) = - lim /EE ajMCjis^k^ds s: 0, 

j—m j=m j—m k—j — m-\-l 

from which we conclude that t \-> Y^^jLm c 3 W * s a mon otonic nonincreasing function. 
As it is bounded below (by zero) it must converge to some constant ^ 0. Since 
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J2T=m c j" W > Ej°l,n+i c i(*) we have P*m > Pm+v Then > for all m e N we have 



OO 



j—m j=m-\-l 

Now consider the case m = 1 and £ = r + 1. Applying limits r — >■ oo to the moments' 
equation, we have 

oo oo pi n oo 

Vcj(t + 1) - Vcj(r) = - lim / V V aj, k Cj(s)ck(s)ds 

* — * ' n— )-oo / ' * 

j=l j=l * /r j = lfe=j 



~ t n oo 

= — lim lim / / / a? kCj(s)ck(s)ds 

r—>oo n-^-oo I ^ — ' ^ — ' ' 
Jt J=lfc=j 



3=1 fc=j 

Suppose there exists an integer p £ N such that c* > 0. Let a p £ (0, c*). Then 

"T + l n 00 







— lim lim / > > a~ kCj(s)ck(s)ds 

t— >oo n^oo I — J — J ' 

lim lim / > a,, 7 (c,(s)) 



/•T+l 2 

^ lim / a p,p(c P (s)) ds > lim a vv o? = a pp a v > 0, 
and this contradiction proves that c* = 0, for all values of p. □ 
9. On the scaling behaviour of solutions 

In this section we begin the study of the scaling behaviour of the solutions of 
(1.1) in the particular case when a,j t k = l,for j, k £ N, in which case the system 
turns into 

00 00 
(9.1) Cj = Y2 c i+kCk-Cj'Y2 Ck > J = 1 ' 2 '--- 

k=l k=l 

This study is strongly motivated by similar studies on the scaling behaviour of 
coagulation-fragmentation equations and other related equations [10]. Most of the 
results in those works are consequences of the application of tools based on the 
Laplace transform. Here we have an entirely new situation since the production 
term (the first one on the r.h.s. of (9.1)) is not of convolution type and hence 
Laplace transform methods are not useful. 

In the first place, we can draw consequences from the differential version of (4.11) 
and (4.12) about the typical time scales for the cluster eating equation. Some of 
these conclusions were already formally obtained in [11] and here we reproduce part 
of their arguments. In that work the authors were led to interesting novel features 
about the evolution of the numbers of clusters of odd and even sizes, which have 
no parallel in the usual Smoluchowski's coagulation-fragmentation equation. These 
were already pointed out in Section 7. Let us define, for t > 0, 

00 00 

K*) : = ]C c -?(*)' ^oddW :=Y^ c 2j-i(t) ■ 
i=i j=i 
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With our choice of the kinetic coefficients, equation (4.11) becomes 

3- 



while equation (4.12) turns into 

(9-3) !>odd = -^odd • 

Equation (9.3) reflects the fact that the number of even size clusters does not affect 
the evolution of the number of odd size clusters. This is intuitively clear since the 
only interaction that changes the net amount of odd sized clusters is the reaction 
between two odd size clusters, which produces an even size one. Solving (9.3) we 
get 

^odd(O) 

^odd (t) = — 77^7, for t > . 

1 + fodd(0)£ 

On the other hand, (9.2) implies that 

-v 2 ^ v < --v 2 , 
2 ' 

and hence 

i/(0) , , v(0) 



1 + v(p)t ^ w " 1 + A 



2 

These results give us a typical time scale of i _1 which was already seen in the 
monodisperse solutions in section 7. From the above computations, 

JO, if ^odd(O) = 0, 
hm t^odd(i) = < 

*"►+«> [1, if ^odd (0)^0, 

liminftj/(t) 1, limsuptz/(t) ^ 2. 

We now turn our attention to the study of self-similar solutions. We call a 
solution c(-) self-similar if there arc functions £, rj, <fi such that, 

(9.4) C ^ ) = ^(i))' te '' eN '*>°- 

We first remark that by virtue of the results obtained in Section 7 with respect to 
the monodisperse solutions (which are a trivial case of self-similar solutions), no 
universal behaviour is to be expected with respect to all the nonnegative solutions 
of (9.1). In other words, for j and t large, the solutions will not have an unique 
asymptotic self-similar behaviour independently of their initial conditions. 

In order to find self-similar solutions we adopt an ansatz suggested in [8] in a 
related situation: find diffcrcntiablc functions A and a such that there is a solution 
in the form 

Cj (t) = A(t)a(ty . 
By plugging this into (9.1) we obtain 

• ■ _ aA 2 

Aa-* + 3A0P a = . 

1 — cr 

This clearly implies that a is a constant in [0, 1), while A must satisfy 

A = --^- 2 A 2 . 

1 — a z 
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By integrating this ODE with initial condition Cj(Q) = Apa^ , for j e N, with 
Aq > 0, a e [0, 1), we obtain the solution 

(9-5) c ^ t) = YTWV J' eN '^ ' 



where (3 := This fits our definition of self-similar solution given by (9.4) 



A a 

I -a 2 
with, 

C(t) = l + pt, »7(t) = l, ^(x)=A a x . 

This confirms again the typical time scale of t^ 1 and, in fact, for each such solution 
and for j G N, 

(9.6) lim te ? (t) = (l-a 2 W -1 , 

t— >+oo 

and 

(9.7) lim = 1 + a. 

t— >-f oo 

Therefore, we can consider (9.5) as a family of self-similar solutions, whose scaling 
behaviour depends only upon the rate a of exponential decreasing of the initial 
condition. For each a £ [0,1) there is one such solution with number of cluster 
density v{t) behaving like in (9.7). At present the answers to the following questions 
about the asymptotic behaviour of the solutions of the cluster eating equation (9.1) 
are still unknown: 

(1) Besides (9.5) are there other strictly positive self-similar solutions? 

(2) Are there other solutions behaving asymptotically like the solutions (9.5) 
in the sense that they satisfy (9.6)? 

(3) If the question above has positive answer and if we know that the total 
number of clusters has the behaviour (9.7) for a known a, can we guar- 
antee that the solution itself behaves like (9.6) or are there other types of 
asymptotic behaviour? 

We hope to return to these questions, as well as to a further exploration of the 
RBK system in the near future. 
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